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I. INTRODUCTION 

Plasma waves are well investigated in condensed matter physics. In the present paper we 
derive the dispersion of plasma waves in the framework of Green functions formalism for T - 0. 
Our purpose is purely pedagogical. We use the diagram technique described in a textbook by 
Lifshitz and Pitaevskii |l|] in order to derive the dielectric polarizability of degenerated electron 
gas. The same result is derived in [2] in a different way. Our analysis starts with the ID case. Then 
we consider the 2D case which is realized in semiconductor inverse layers. For a review of the 
properties of two dimensional electron gas see, for example, Ref . |g. 

II. MODEL 

The dielectric formalism gives us the equation for the dielectric polarizability 

e / (Q,q) = l-y(Q,q)n(Q / q). (1) 

In that expression V(C1, q) means photon propagator, which in non-relativistic approximation 
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is the Fourier transformation of the Coulomb potential 

V*»(Q,q)«— , ^q)^. (2) 

The polarization operator II is the functional derivative of the electron density with respect to 
small perturbations of the electric potential. The diagram rules described in Ref. Ill], section 13 
give for the simple electron loop 

d D+1 P 

where G a p(P) are the Green functions of degenerated electron gas 



n(Q,q) = (-O(-l) JiG^DiGf^P + Q)-^^, (31 



We follow standard notations: r\ p - e p - ii, /i is the Fermi energy, e p = p 2 /2m is the electron kinetic 
energy and P = (<u, p), Q - (Q, q). The dispersion of longitudinal plasma waves Cl(q) is the 
solution of 

£/(Q,q) = 0. (5) 

The first step is calculating the polarization operator. Using 

G(P)G(P + Q) = . ^ 1 = (6) 

co - T] p (l -i0) co + Q- rj p+q (l - zO) 
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1]p ^Ip+q + ^ 



cu - j] p (l -z'O) cu + Q - T] p+(? (1 - z'O) 



the cu integration gives the electron filling factors 



CO - 7] p (l - Z'O) 

In such a way we obtain 

d D p C doo 



d°P "p - "p+q 



= Il+I 2 , (8) 



(27l) D 7]p - T] p+q + Q + z'O 

where the evanescent imaginary part of the frequency corresponds to the Landau rule for adiabatic 
appearance of the perturbation, see Ref. [2|] sees. 29, 40, Eqs. (29.6), (40.15). Here for brevity we 
use the notations 

f d °P n P _ f d D p n p 

7l(Q ' q)_2 J (27i)Oz ] p-z ] p +(? + Q- 2 J (271)° v±_^t^ (9) 
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j (O \ = ? f d V ~" p+t? =2 f d V ~ np+? 
2 ^ ' q) J (Inpnp-rim + Q J (2n)Vf_(m£ + Q 
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We change the variables p' — p + q. Then dp — dp' , n v+ q — rip> and 
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(11) 



2m 2m 

The general formula © for the response function of the electron gas is given in many textbooks 
on statistical and solid state physics. 

R 

In the next section we will analyze the simple case of D = 1 which is discussed in [3]. 



III. 1-DIMENSIONAL CASE 



For the case of D = 1 we assume q > 0, and 
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The velocity on Fermi level is v F - p F /m, and 
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Owing to the symmetry of the polarization operator 

n(Q) = n(-Q), 
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we obtain 
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In hydrodynamic approximation (long waves and phase velocity much larger than Fermi velocity) 

(17) 
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we get 
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Let us introduce D-dimensional density of states 

T7 TT D / 2 

" (DU2 (2^' W= (D72)! pD ' ^(PF ) = 2 PF for D = 1, 
where Vp is the volume of the Fermi sphere. In our case 

2n n 
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Then 

q 2 n (W) 

n = - — (22) 

Calculating the Fourier transformation, it can be shown ] 3, 4 , ^] that the photon propagator in 
long wavelength approximation (q — > 0) obeys 

V( 1D >(Q /i? )«2 e 2 ln|iJ, (23) 

where a> , a — » is a regularization parameter. More precisely 

V^ D \Q, q ) = e 2 f ^. (24) 

This integral is not convergent - we have to regularize it introducing a small diameter of the wire. 
We give two schemes: 

1) V( w \n,q) = e 2 ^ j=j = -^(No(iaq) + N (-iaq))*2e 2 ]n(— ), (25) 

2) y( 1D )(Q, q) = le 1 f°° cos ^ dx = - 2e 2 Ci(qa) « -2e 2 (y + Infafl)) = 2e 2 In (— 

Jfl x \eYqa 

Here y is Euler's constant. As it is written above, both schemes yield 

V «2e 2 ln| — 



(26) 



In order to find the dispersion we have to solve the equation (O 
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As a result, for dispersion of plasma waves we obtain 



le 2 2 n (lD) / 1 

VU= — 2__]n _ =i. (27) 



Q 2 (?) = ^— ^ ^ U 2 , (28) 
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in agreement with the derivation in |3fl. 

In the next section we will discuss the case of D - 2 which is for first time was discussed in (6fl. 

IV. 2-DIMENSIONAL CASE 

For the case of D = 2 we introduce polar variables 

p = p(cos 6, sin 9), q = q(l, 0), 
d 2 p = pdpdO, p • q = pq cos 6, 



and 
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The velocity on Fermi level is vp = pp /m, and 



nr 



I 1 (Q,q) = — 
uq L 



Q - — sign(Q - — ) -i /( Q - — ) - q 2 vl 



2m' 



1 
2m 



m 



h(Q,q) = — 
nq L 



-Q - — + sign(Q + — ) -i /(q + — 1 -n2 T r- 



2m' 



2m 



Owing to the symmetry of the polarization operator 



n(Q) = n(-Q), 



we obtain for big enough frequency 
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Working in hydrodynamic approximation ((T7)) , we introduce the notation 
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In these terms we have for the polarization operator 
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we arrive at 
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For the 2D case the electron density 
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Then 
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Again for the dispersion we solve the equation (|5]) 



vn = — n — — = i. 
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As a result for dispersion of 2D plasma waves, we get 
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Another important case is the static one (Q = 0) for 2D electron gas where for q < 2pp we have 

Re(n) = -s B]. 

In the next section we will discuss the D = 3 case which is described in the monograph by 
Abrikosov, Gor'kov and Dzyaloshinsky Jz]. 



V. 3-DIMENSIONAL CASE 



Analogously for D = 3 we introduce spherical coordinates 



p = p(sin cos cp, sin sin <p, cos 6), q = ^(0, 0, 1), 
d 3 p = p 2 sin dp d6 d<p, p • q = p q cos 6. 
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Then (we are not interested in the imaginary part) 
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To simplify the notations we introduce a = 2raQ — q 2 , b = 2q, so 
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and 
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Finally we obtain 
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The Fourier transformation of the static polarization operator gives the we 
physics of magnetism RKKY interaction. Technical details are given in Ref . {]. 
Introducing the notations 

m(Q 2 - q 2 vl) Ci + qv F , 

4 3 3 

Q e = V47in(3D) e 2/ m , „(3D) = 2 |!Z| = £F / 
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the longitudinal polarizability takes the form J2] 
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As we did before, we use the hydrodynamical approximation ((17)) , and 
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That is why equation (|45l ) transforms into 
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4mp 3 q 2 / 1 1 \ P 3 ^ 2 

+ 3ti 2 \(2mQ-^ 2 ) 3 ~ (2mQ + q 2 f J ~ 3n 2 mD 2 ' 

Taking into account the density of states we can write the polarization operator in another form 

(51) 

This hydrodynamic asymptotic is valid for arbitrary dimensions. 
Turning back in (0 

^yW 51 (52) 

we finally obtain the well-known result for bulk plasma oscillations 

47ze 2 «( 3D ) 

Q 2 (o) = = const. (53) 

m 
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